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Abstract 

Inclusive jet cross sections in photoproduction for events containing a D* meson 
have been measured with the ZEUS detector at HERA using an integrated lumi- 
nosity of 78.6 ptT 1 . The events were required to have a virtuality of the incoming 
photon, Q 2 , of less than 1 GeV 2 , and a photon-proton centre-of-mass energy in 
the range 130 < W w < 280 GeV. The measurements are compared with next- 
to-leading-order (NLO) QCD calculations. Good agreement is found with the 
NLO calculations over most of the measured kinematic region. Requiring a sec- 
ond jet in the event allowed a more detailed comparison with QCD calculations. 
The measured dijet cross sections are also compared to Monte Carlo (MC) mod- 
els which incorporate leading-order matrix elements followed by parton showers 
and hadronisation. The NLO QCD predictions are in general agreement with 
the data although differences have been isolated to regions where contributions 
from higher orders are expected to be significant. The MC models give a better 
description than the NLO predictions of the shape of the measured cross sections. 
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1 Introduction 



Charm and/or jet production in ep collisions should be accurately calculable in pertur- 
bative Quantum Chromodynamics (pQCD) since the mass of the heavy quark, tuq, and 
the transverse energy of the jet, E^ , provide hard scales. In photoproduction, where 
a quasi-real photon, emitted from the incoming lepton, collides with a parton from the 
incoming proton, such events can be classified into two types of process in leading-order 
(LO) QCD. In direct processes, the photon couples as a point-like object in the hard 
scatter. In resolved processes, the photon acts as a source of incoming partons with only 
a fraction of its momentum participating in the hard scatter. 

Measurements of the D* photoproduction cross section [1] as functions of the transverse 
momentum, p®* , and the pseudorapidity, t] D * , show that the predictions from next-to- 
leading-order (NLO) QCD are too low for p®* > 3 GeV and t] D * > 0. Part of this deficit 
may be due to hadronisation effects. The predictions for jet production accompanied by a 
D* meson should have smaller uncertainties from these hadronisation effects. Furthermore 
jets can be measured in a wider pseudorapidity range than D* mesons due to the larger 
acceptance of the calorimeter compared to the central tracker. 

A dijet sample of D* photoproduction can also be used to study higher-order QCD topolo- 
gies [1,2]. In the present paper, previously unmeasured correlations between the two jets 
of highest transverse energy, namely the difference in azimuthal angle, A(/$, and the 
squared transverse momentum of the dijet system, (p 1 ^) 2 , which are particularly sensitive 
to higher-order topologies, are presented. For the LO 2 — > 2 process, the two jets are 
produced back-to-back with A^ = ir and (p?) 2 = 0. Large deviations from these values 
may come from higher-order QCD effects. The accuracy of the theoretical description of 
these effects is tested. 

Calculations performed to NLO in QCD are available with two different treatments for 
charm. In the fixed-order, or "massive", scheme [3], u, d and s are the only active flavours 
in the structure functions of the proton and photon; charm and beauty are produced only 
in the hard scatter. This scheme is expected to work well in regions where the transverse 
momentum of the outgoing c quark is of the order of the quark mass. At higher transverse 
momenta, the resummed or "massless" scheme [4,5] should be applicable. In this scheme, 
charm and beauty are regarded as active flavours (massless partons) in the structure 
functions of the proton and photon and are fragmented from massless partons into massive 
hadrons after the hard process. 

In this paper, photoproduction of charm is studied by tagging a D* ± meson and recon- 
structing at least one jet in the final state. The measurement is performed in the following 
kinematic region: photon virtuality, Q 2 < 1 GeV 2 ; photon-proton centre-of-mass energy, 
130 < W 1P < 280 GeV; p£* > 3 GeV; \r] D *\ < 1.5; jet transverse energy, > 6 GeV; and 
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jet pseudorapidity, —1.5 < 77^ < 2.4. Differential cross sections as a function of E^ and 
rf ct have been measured. Jets are divided into two categories: jets of the first category 
are associated with the D* meson (D*-tagged jet), while jets of the second category are 
not matched to a D* meson (untagged jet). The inclusive, L>*-tagged and untagged jet 
cross sections are compared to the massive NLO QCD predictions. A comparison to the 
massless calculation is only available for the untagged jet cross sections [6]. 

A sub-sample having at least two jets with > 7GeV and E^ t2 > 6 GeV is used to 
measure the correlations between the two highest Ej, jets: the fraction of the photon 
momentum participating in dijet production, x° bs [7]; A0jj; (p^) 2 ; and the dijet invariant 
mass, M J j. Differential dijet cross sections as a function of these variables have been mea- 
sured for the direct-enriched (x° bs > 0.75) and resolved-enriched (x° bs < 0.75) kinematic 
regions and compared to massive NLO QCD predictions and Monte Carlo (MC) models. 

2 Experimental set-up 

The analysis was performed with data taken from 1998 to 2000, when HERA collided 
electrons or positrons 1 with energy E e = 27.5 GeV with protons of energy E p — 920 GeV 
resulting in a centre-of-mass energy of 318 GeV. The results are based on an integrated 
luminosity of 78.6 pb" 1 of ep collision data taken by the ZEUS detector. A detailed 
description of the ZEUS detector can be found elsewhere [8]. A brief outline of the 
components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [9], which oper- 
ates in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD 
consists of 72 cylindrical drift chamber layers, organised in 9 superlayers covering the 
polar- angle 2 region 15° < 9 < 164°. The transverse-momentum resolution for full-length 
tracks is a(p T )/p T = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [10] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a(E)/E = 0.18 /y/E for electrons 
and a{E)/E = 0.35/ \[E for hadrons (E in GeV). 

1 Hereafter, both electrons and positrons are referred to as electrons, unless explicitly stated otherwise. 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the center of HERA. The coordinate origin is at the nominal interaction point. 
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The luminosity was measured from the rate of the bremsstrahlung process ep — > e7p, 
where the photon was measured in a lead-scintillator calorimeter [11] placed in the HERA 
tunnel at Z = —107 m. 

3 Event reconstruction and selection 

A three-level trigger system was used to select events online [8, 12]. At the first- and 
second-level triggers, general characteristics of photoproduction events were required and 
background due to beam-gas interactions rejected. At the third level, a D* candidate was 
reconstructed. 

In the offline analysis, the hadronic final state and jets were reconstructed using a combi- 
nation of track and calorimeter information that optimises the resolution of reconstructed 
kinematic variables [13]. The selected tracks and calorimeter clusters are referred to as 
Energy Flow Objects (EFOs). To select photoproduction events, the following criteria 
were used: 

• the event vertex was required to be within 50 cm of the nominal vertex position in the 
longitudinal direction; 

• deep inelastic scattering (DIS) events with a scattered electron candidate in the CAL 
were removed [14]. To keep the events where a pion was misidentified as a scattered 
electron, events where y e = 1 — ^-(1 — cos9 e /) > 0.7 were retained; E e > and 9 e > are 
the energy and polar angle, respectively, of the scattered electron candidate; 

• the requirement 130 < Wjb < 280 GeV was imposed, where Wjb = a/ ^E p E e y jb and 
Ujb is the estimator of the inelasticity, y, measured from the EFOs according to the 
Jacquet-Blondel method [15]. Here, Wjb was corrected, using MC simulation, for the 
energy losses of EFOs in inactive material in front of the CAL. The upper cut removed 
DIS events where the scattered electron was not identified and which, therefore, have 
a value of yjB close to 1. The lower cut removed proton beam-gas events which have 
a low value of yjB- 

The cuts on y e and Wjb restricted the range of the virtuality of the exchanged photon to 
Q 2 less than about 1 GeV 2 , with a median value of about 3 x 10 -4 GeV 2 . 

3.1 D* reconstruction 

The D* mesons were identified using the decay channel D* + — > D°ir^ with the subsequent 
decay D° — > K~ir + and the corresponding antiparticle decay, where 7r+ refers to a low- 
momentum ("slow") pion accompanying the D°. 
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Charged tracks measured by the CTD and assigned to the primary event vertex were 
selected. The transverse momentum was required to be greater than 0.12 GeV. Each 
track was required to reach at least the third superlayer of the CTD. These restrictions 
ensured that the track acceptance and momentum resolution were high. Tracks in the 
CTD with opposite charges and transverse momenta px > 0.4 GeV were combined in pairs 
to form D° candidates. The tracks were alternately assigned the masses of a kaon and a 
pion and the invariant mass of the pair, M^, was evaluated. Each additional track, with 
charge opposite to that of the kaon track, was assigned the pion mass and combined with 
the D°-meson candidate to form a D* candidate. 

The signal regions for the reconstructed masses, M(D°) and AM = (Mkwit 3 — Mk-k), were 
1.80 < M(D°) < 1.92 GeV and 0.143 < AM < 0.148 GeV, respectively. For background 
determination, D° candidates with wrong-sign combinations, in which both tracks forming 
the D° candidates have the same charge and the third track has the opposite charge, were 
also retained. The same kinematic restrictions were applied as for those D° candidates 
with correct-charge combinations. The normalisation factor of the wrong-charge sample 
was determined as the ratio of events with correct-charge combinations to wrong-charge 
combinations in the region 0.15 < AM < 0.17 GeV. 

The kinematic region for D* candidates was p®* > 3 GeV and \rj D *\ < 1.5. Figure Q 
shows AM for the selection of a D* meson with a jet (see Section f3.2|) . The fit to the 
distribution has the form 

F = pi ■ exp (-0.5 • x 1+ ) + p 4 ■ (AM - m w ) P5 , 

where x = |(AM— P2)/ps\, Pi~P5 are free parameters and m» is the pion mass. The "mod- 
ified" Gaussian described both data and MC distributions well. The fit gives a peak at 
145.467 ±0.015(stat.) MeV to be compared to the PDG value of 145.421 ± 0.010 MeV [16]. 
The difference is due to systematic effects which are too small to be relevant for this anal- 
ysis. The fitted width of 0.61 ± 0.02 MeV is consistent with the experimental resolution. 

The number of D* mesons was determined from candidates reconstructed in both signal 
regions and after the subtraction of the background estimated from the wrong-charge 
sample; this gave 4891 ± 113 D* mesons. This procedure was used throughout the paper 
with the number of D* mesons obtained from the fit used as a systematic check. 

3.2 Jet reconstruction 

Jets were reconstructed using EFOs as input to the the kj< cluster algorithm [17] in its 
longitudinally invariant inclusive mode [18]. The transverse energy of the jet was corrected 
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for energy losses in inactive material in front of the CAL, where the correction factors 
were determined in bins of E^ and 7f et from MC simulation. These corrections were 
between 5% and 12%. 

For the inclusive jet cross sections, jets with > 6 GeV and —1.5 < rf ct < 2.4 were 
selected, and events containing at least one such jet were used for further analysis. For the 
dijet analysis, events were required to have at least two jets with —1.5 < rf et < 2.4 and 
E^ > 6 GeV, while the highest E^ jet was required in addition to satisfy E^ 1 > 7 GeV. 
The asymmetric jet transverse-energy cut assures that the NLO calculation is not infrared 
sensitive [19]. After the dijet selection, 1692 ± 70 D* mesons remained. 

Cross sections are given separately for Z)*-tagged and untagged jets. _D*-tagged jets are 
defined as jets in which a D* (in the kinematic region defined in Section I3~T]) was clustered 
into the jet at the hadron level by the kx algorithm. All the other jets are called untagged 
jets. These two classes of jets can be distinguished experimentally by cutting on the 
distance AR(D*, jet) = (77-°* — ?7j et ) 2 — (0 D * — ^ ct ) 2 between the D* and the jet, where 
<p D * and 0-> et are the azimuthal angles of the D* meson and jet, respectively. Figure 121 
shows the distance of the D* to all jets in the event for data and MC. The peak at 
AR(D*,jet) = is due to _D*-tagged jets as indicated by the Ai?(D*,jet) distribution for 
-D*-tagged jets for the MC hadron-level predictions also shown in Fig. 121 The broad peak 
at AR(D*, jet) ~ 3 is due to the untagged jets. A cut at AR(D*, jet) < 0.6 was used to 
distinguish experimentally tagged and untagged jets. For the inclusive jet sample, 83% 
of D* mesons were matched to a jet, and for the dijet sub-sample 94% of D* mesons were 
matched to a jet. 

4 Monte Carlo models 

The MC programmes Herwig 6.301 [20,21] and Pythia 6.156 [22], which implement 
LO matrix elements followed by parton showers and hadronisation, were used to model 
the final state. The Herwig and Pythia generators differ in the details of the im- 
plementation of the leading-logarithmic parton-shower models. They also use different 
hadronisation models: Herwig uses the cluster [23] model and Pythia uses the Lund 
string [24] model. Direct and resolved events were generated separately and in proportion 
to the cross sections predicted by the MC programme. The relative fraction of charm and 
beauty events was also generated in proportion to the cross sections predicted by the MC 
programme. Events were generated using CTEQ5L [25] and GRV-G LO parton density 
functions (PDF) for the proton and the photon, respectively. The c-quark and 6-quark 
masses were set to m c = 1.5 GeV and mj = 4.75 GeV, respectively. 

The MC programmes were used both to compare with the dijet cross sections, which are 
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particularly sensitive to the parton-shower models and for calculation of the acceptance 
and effects of detector response (see Section EJ). For all generated events, the ZEUS 
detector response was simulated in detail using a programme based on Geant 3.13 [26]. 

5 NLO QCD calculations 

There are two NLO QCD calculations available to calculate jet cross sections in charm 
photoproduction: the massive calculation by Frixione et al. (FMNR) [3] and the massless 
calculation by Heinrich and Kniehl [6]. 

5.1 Massive calculation 

In the massive calculation, the PDF sets used were CTEQ5M1 [25] for the proton and 
AFG-HO [27] for the photon. The renormalisation scale, /xr, and factorisation scale, /xp, 
were set to /x = /x# : = jip = = \J ((Pt) 2 } + m c> where ((pt) 2 } is the average squared 
transverse momentum of the two charm quarks and m c = 1.5 GeV. The fragmentation 
of the charm quark into a D* meson was described by rescaling the c-quark momentum 
using the Peterson fragmentation function [28] with e = 0.035 ±0.002 which is taken from 
an NLO fit to ARGUS data [29]. The fraction of charm quarks hadronising into a D* 
meson was set to 0.235 [30]. The kr algorithm was applied to the outgoing partons in the 
final state of the NLO programme. 

The dependence of the NLO prediction on different photon PDFs, fin, \ip and m c was 
evaluated by repeating the calculation using different sets of parameters. The upper 
(lower) bound of the NLO QCD prediction was estimated by setting fiR = rriT/2 and 
m c = 1.3 GeV (/xr = 2m*r and m c = 1.7 GeV). 

An NLO prediction of D* production for beauty is not available so this contribution was 
estimated using a combination of the B hadron cross section at NLO and B decays in 
Pythia. The pt distributions of the two stable B hadrons produced in the Pythia 
MC programme were reweighted to the distribution in the NLO calculation. In the 
NLO calculation, the 6-quark mass, m&, was set to 4.75 GeV, /x = = a/ ((Pt) 2 ) + m b 
and e = 0.0035 [31]. The branching b — > D* was set to the value measured by the 
OPAL Collaboration [32]. The upper (lower) bound of the NLO QCD prediction was 
estimated by setting /i R = m T /2 and m h = 4.5 GeV (jar = 2m T and m c = 5.0 GeV). 
The contribution from beauty production for the inclusive jet distribution, as predicted 
by NLO+Pythia, is about 2% at low and increases to 8% at high E^ . For each 
cross section, this beauty contribution and its uncertainty were added linearly to the 
corresponding massive D* prediction from charm quarks. 
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5.2 Massless calculation 



In the massless calculation, AFG04 [33] for the photon PDF and MRST03 [34] for the 
proton PDF were used. The number of flavours was set to five. The D* fragmentation 
function and fraction of beauty and charm hadronising into a D* meson are derived [5,6] 
from a fit to data from LEP; the function is assumed to be applicable to HERA as it is 
derived using the factorisation theorem in QCD. The central prediction uses lir = m' T = 
V (Pt*) 2 + m c> where m c = 1.5 GeV. The fragmentation factorisation [6] scale, Mp, and 
fip were set to Mp — fip = 2m' T . 

The uncertainty was estimated by changing the scales to lir = m' T /2 and \xp = Mp = Am' T 
for the upper bound, and lir = 2m' T and /ip = Mp = m' T for the lower bound. The 
photon PDF, GRV-HO, was used, and the uncertainty was found to be less than half that 
from the variation in scale [6]. The difference between MRST01 and MRST03 proton 
PDFs was found to be negligible for the distributions considered. The size of the beauty 
contribution was estimated by suppressing the final-state fragmentation of a b quark to 
a D* meson. This reduces the cross-section prediction by about 3% at low and 15% 
at high E^. Due to theoretical limitations, the predictions for the massless scheme can 
only be calculated for the untagged-jet distributions. 



5.3 Hadronisation correction 

As the NLO calculations produce final-state partons, the effects of hadronisation were 
considered when comparing the predictions with the data. The NLO QCD jet predictions 
were corrected using a bin-by-bin procedure according to da = <icx NLO - Chad, where <icr NLO 
is the cross section for parton jets in the final state of the NLO calculation. The hadronisa- 
tion correction factor was defined as the ratio of the jet cross sections after and before the 
hadronisation process, Chad = da^ rons /da^c tons . Here, parton-level cross sections were 
obtained using partons after the initial- and final-state showering of the MC simulations 
described in Section |3J Distributions at the parton level in the MC programmes were 
checked to be similar to those calculated using the NLO programme, assuring the validity 
of using a bin-by-bin correction. The value of Chad was taken as the mean of the ratios 
obtained using the Herwig and Pythia predictions. The uncertainty on this value was 
estimated as half the difference between the values obtained using the two models. These 
uncertainties were added in quadrature to the other uncertainties of the NLO calculations. 
The values of Chad applied to the NLO predictions are given in Tables fTHTTl 
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6 Data correction and systematic uncertainties 



The data were corrected, using the MC models described in Section 0J for the detector 
acceptance and the selection efficiencies to obtain differential cross sections for the process 
ep — > el + D* + jet(s) + X. The definition of the cross section includes events with a 
D* containing primary c quarks or those from 6-quark decays. The data were initially 
compared to the MC simulation in shape and found generally to agree well for all the 
kinematic quantities. Since Herwig gives a better overall description of the data than 
Pythia, it was chosen as the primary MC model to correct the data. The cross section 
for a given observable Y was determined using 

da _ N 

dY ~ A-C-B-AY' 

where N is either the number of jets for the inclusive jet cross sections or the number of 
D* mesons for the dijet cross sections in a bin of size Ay. The acceptance, A, takes into 
account migrations and efficiencies for that bin, and £ is the integrated luminosity. The 
product, B, of the appropriate branching ratios for the D* and D° was set to (2.57 ± 
0.06)% [16]. 

The systematic uncertainties of the measured cross sections were determined by changing 
the selection cuts or the analysis procedure in turn and repeating the extraction of the 
cross sections. The uncertainties are described in more detail elsewhere [35]. The following 
systematic studies have been carried out (the resulting uncertainty on the total cross 
section is given in parentheses): 

• varying the values of the selection cuts by the experimental resolutions in the corre- 
sponding quantity (1^3%); 

• varying the efficiencies of the CAL first-level trigger (+4.1%); 

• the acceptance was recalculated by re-weighting the prediction from the Herwig 
MC simulation in pj,* / to reproduce the distribution of this variable in the data 
(+5.0%); 

• an additional contribution to the uncertainty from the modelling of the hadronisation 
process was estimated by using Pythia instead of Herwig (+0.4%); 

• the effect of the uncertainty of the beauty cross section on the acceptance correction 
was taken into account by increasing the beauty contribution by a factor two (< ±1%); 

• varying the procedure to extract the D* signal {-2.9%); 

• varying by ±3% the jet energy scale in the CAL (tl'1,%)- 
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All systematic uncertainties were added in quadrature to obtain the total systematic 
uncertainty, except for the jet energy-scale uncertainty, as this has a large correlation 
between bins, and is shown separately in all figures. In most bins of the differential 
cross sections, the total systematic uncertainty is comparable to the statistical errors. In 
addition, an overall normalisation uncertainty of 2% from the luminosity determination 
is included in neither the figures nor the tables. 

7 Inclusive jet cross sections 

Inclusive jet cross sections with a D* in the final state have been measured as a function 
of .Ej^ and 7f ct in the following kinematic region: 

• Q 2 < 1 GeV 2 and 130 < W 1V < 280 GeV; 

• v¥ > 3 GeV and |^ D *| < 1.5; 

• E^ > 6 GeV and -1.5 < rf ct < 2.4; 

• a -D*-jet match is required only where explicitly stated. 

At the hadron level, the D* meson is used as input into the jet finder therefore allowing 
a £>*-tagged jet to be identified unambiguously. 

The cross-sections da/dE^ in bins of rf ct are shown in Fig. Eland given in Tabled for all 
jets and in Fig. H]and Table El for the whole rf ct range for D*-tagged and untagged jets. 
The distributions tend to fall less steeply with increasing rf et . The cross sections da/drf ct 
in different regions of E^ are shown in Fig. El and given in Table 121 for all jets and in 
Fig. El and Tables ?? and ?? for _D*-tagged and untagged jets. Due to the requirement 
\rj D * | < 1.5, the -D*-tagged jet is centred around rf et = and falls off rapidly at large rf ct . 
The advantage of reconstructing jets is observed in the untagged-jet distribution where a 
significant cross section is measured up to rf et = 2.4. This is due to the larger acceptance 
of the CAL compared to the CTD. 

The massive calculation is compared to all measured cross sections whereas the massless 
calculation is compared only to the untagged-jet distributions. The normalisation of the 
data for all distributions is well described by the upper limit of both NLO QCD predic- 
tions. The shape of the data is well described by the NLO QCD predictions. The inclusion 
of hadronisation corrections which shift the distributions towards forward rf et improves 
the description of the data. Some difference in shape is observed for the upper bound of 
the massless prediction compared to the untagged-jet distributions for E^ > 9 GeV (see 
Fig. ED. 

Measurements of dajdE^ and da/drf ct in bins of p®* are given in Tables El and El and 
shown in Fig. [7| compared to the massive NLO prediction. The NLO prediction gives a 
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poor description of the normalisation of the data at lowest p® ■ However, the normalisa- 
tion of the NLO prediction agrees with the data in the two regions of higher p^* ■ In all 
regions, the shape of the data is reasonably well described by the NLO prediction. Similar 
conclusions on the normalisation were seen for inclusive D* measurements [1]. However, 
the difference in shape observed as a function of r] D * in the inclusive measurement is not 
seen here as a function of rf ct . The 6-quark contribution is largest in the region of low 

vW- 

In order to be sensitive to higher-order effects, and to distinguish between direct-enriched 
and resolved-enriched regions, the variable x° bs (.D*, jet) was constructed [6], in an analo- 
gous way to the 'traditional' x° bs [7]. Using the D* meson and the untagged jet of highest 
Ej?, the quantity x° bs (D*, jet) is given by: 

„D*„-ri D * , 7-i(untagged jet) „(untagged jet) 

= ^ +E - 2yEt . (i) 

This variable has the advantage of being calculable in the massless scheme. In addition 
it takes advantage of increased statistics by requiring only one jet of high E^ . In Fig. ED 
the measured cross-section da/dx^ 3S (D*, jet), given in Tabled is compared to both the 
massless and massive predictions. The upper bound of the massive prediction gives a 
good description of the data; the description of the massless prediction is somewhat 
worse. The Herwig MC model gives a poor description whilst Pythia gives a reasonable 
description of the shape of the data distribution. Both MC programmes underestimate 
the normalisation of the data. 



8 Dijet cross sections 

Dijet correlations are particularly sensitive to higher-order effects and therefore suitable 
to test the limitations of fixed-order perturbative QCD calculations. Events containing a 
D* meson were required to have at least two jets with E^ n > 7 GeV, E^ t2 > 6 GeV and 
— 1.5 < rf et1 ' 2 < 2.4. The Q 2 , W w , p®* and r] D * requirements were the same as for the 
inclusive jet cross section. 

The dijet variables measured were reconstructed from the two highest E 3 ^ jets as: 



pjetl ^ti R jet2 -jfo 
nbs _ T ' T (n\ 

i ~ 2yE e ' {2) 

A0 jj = |0 jctl -<jj ct2 \ , (3) 
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(4) 2 = (pT 1 + 2 ) 2 + (pT + pf 2 ) 2 , (4) 



M jj = ^2E$ n E% t2 [cosh(r/j etl - r/j ct2 ) - cos(<^ ctl - & ct2 )} . (5) 

Table El gives, and Fig. shows, the dijet cross section as a function of x° bs , which is 
reasonably well described by the massive calculation. In Tables CTTT1 and Figs. H3>d, the 
cross sections as a function of A(f$, {p^) 2 and M ]i are also shown. For A(f$ there is agree- 
ment between data and the NLO prediction at large angular separation, but at smaller 
values the NLO prediction underestimates the data. This is correlated with the 
agreement and disagreement at low and high (p^) 2 values, respectively. The distribution 
in dijet invariant mass is described well by the upper limit of the NLO prediction, as was 
the case for the inclusive jet cross sections in Section [7| 

Cross sections as a function of M 1 *, A^ and (p^) 2 are shown in Figs. ll0fJT2l and given in 
TableslMTTlseparatelv for direct-enriched (x° bs > 0.75) and resolved-enriched (x° bs < 0.75) 
samples. The data are compared to massive NLO QCD predictions and expectations from 
MC models. 

The cross-section da/dM^ in Fig. ^] is described well by the NLO prediction and both 
MC models, Herwig and Pythia, for both regions in x° bs . 

The cross-sections da/dAft* (see Fig. ITTj) and da/d(p^) 2 (see Fig. IT2*j) are reasonably 
well reproduced by the NLO prediction for x° bs > 0.75 although the data exhibit a 
somewhat harder distribution. For x° bs < 0.75, the data exhibit a harder spectrum than 
for x° bs > 0.75. The NLO prediction of the cross section for x° bs < 0.75 has a significantly 
softer distribution compared to the data, both as a function of A(j$ and (pr) 2 - The 
low-x° bs region is more sensitive to higher-order topologies not present in the massive 
NLO calculation. The predictions from the Pythia MC reproduce neither the shape 
nor the normalisation of the data for low and high x° bs . However, the predictions from 
the Herwig MC give an excellent description of the shapes of all distributions, although 
the normalisation is underestimated by a factor of 2.5. The fact that a MC programme 
incorporating parton showers can successfully describe the data whereas the NLO QCD 
prediction cannot indicates that the QCD calculation requires higher orders. Matching of 
parton showers with NLO calculations such as in the MC@NLO programme [36], which 
is not currently available for the processes studied here, should improve the description 
of the data. 
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9 Conclusions 



Differential inclusive jet cross sections for events containing a D* meson have been mea- 
sured with the ZEUS detector in the kinematic region Q 2 < 1 GeV 2 , 130 < W 1P < 
280 GeV, > 3 GeV, \r] D *\ < 1.5, > 6 GeV and -1.5 < rf ct < 2.4. The measure- 
ments are compared to NLO QCD predictions in the massive and massless schemes. With 
the addition of hadronisation corrections, the upper limit of both theoretical predictions 
show similar trends and reasonable agreement with all measured cross sections. Dijet 
correlation cross-sections da/dM ]i and da/dx^ are well described by the massive NLO 
QCD prediction, although again the data tends to agree better with the upper bound 
of the NLO calculation. In contrast, the cross-sections da/dA^ and da/d(p!^) 2 show a 
large deviation from the massive NLO QCD prediction at low A(j$ and high (py*) 2 . This 
discrepancy is enhanced for the resolved-enriched (x° bs < 0.75) sample. These regions 
are expected to be particularly sensitive to higher-order effects. The Herwig MC model 
which incorporates leading-order matrix elements followed by parton showers and hadro- 
nisation describes the shape of the measurements well. This indicates that for the precise 
description of charm dijet photoproduction, higher-order calculations or the implementa- 
tion of additional parton showers in current NLO calculations are needed. 
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range 




range 


da j cLE^ 


± stat. 


± syst. 


± E-scale (nb/GeV) 


-1.5 


2.4 


6.0 


9.0 


1.583 


±0.039 


+0.145 
-0.079 


+0.072 
-0.060 






9.0 


13.0 


0.367 


±0.018 


+0.038 
-0.017 


+0.025 
-0.023 






13.0 


18.0 


0.0775 


±0.0092 


+0.0147 
-0.0073 


+0.0097 
-0.0074 






18.0 


25.0 


0.0269 


±0.0046 


+0.0053 
-0.0055 


+0.0032 
-0.0026 


-1.5 


-0.5 


6.0 


9.0 


0.321 


±0.017 


+0.032 
-0.010 


+0.025 
-0.018 






9.0 


13.0 


0.0408 


±0.0062 


+0.0053 
-0.0039 


+0.0044 
-0.0044 






13.0 


18.0 


0.0088 


±0.0033 


+0.0013 
-0.0055 


+0.0020 
-0.0022 






18.0 


25.0 


0.00 


±0.00 


+0.00 
0.00 


+0.00 
0.00 


-0.5 


0.5 


6.0 


9.0 


0.669 


±0.023 


+0.064 
-0.033 


+0.027 
-0.023 






9.0 


13.0 


0.159 


±0.011 


+0.017 
-0.010 


+0.011 
-0.010 






13.0 


18.0 


0.0315 


±0.0050 


+0.0056 
-0.0030 


+0.0042 
-0.0037 






18.0 


25.0 


0.0089 


±0.0025 


+0.0032 
-0.0012 


+0.0012 
-0.0010 


0.5 


1.5 


6.0 


9.0 


0.446 


±0.022 


+0.048 
-0.024 


+0.016 
-0.013 






9.0 


13.0 


0.122 


±0.011 


+0.015 
-0.010 


+0.007 
-0.006 






13.0 


18.0 


0.0306 


±0.0063 


+0.0110 
-0.0032 


+0.0034 
-0.0020 






18.0 


25.0 


0.0126 


±0.0033 


+0.0032 
-0.0047 


+0.0013 
-0.0014 


1.5 


2.4 


6.0 


9.0 


0.151 


±0.017 


+0.015 
-0.026 


+0.004 
-0.006 






9.0 


13.0 


0.0462 


±0.0072 


+0.0075 
-0.0058 


+0.0024 
-0.0024 






13.0 


18.0 


0.0095 


±0.0041 


+0.0055 
-0.0047 


+0.0013 
-0.0008 






18.0 


25.0 


0.0053 


±0.0019 


+0.0013 
-0.0015 


+0.0007 
0.0002 



Table 1: The cross section da/dE^ for events containing at least one D* meson 
in different regions of rf et . The statistical (stat.), systematic (syst.) and energy 
scale (E-scale) uncertainties are shown separately. 
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Eyi range 




range 


da/drf et 


± stat. 


± syst. 


± E-scale (nb) 


> 6 


-1.5 


-1.0 


0.548 


±0.069 


+0.086 
-0.156 


+0.065 
-0.053 




-1.0 


-0.5 


1.745 


±0.094 


+0.210 
-0.059 


+0.135 
-0.100 




-0.5 


0.0 


2.79 


±0.12 


+0.26 
-0.18 


+0.16 
-0.14 




0.0 


0.5 


2.90 


±0.12 


+0.31 
-0.11 


+0.13 
-0.12 




0.5 


1.0 


2.30 


±0.12 


+0.28 
-0.10 


+0.11 
-0.09 




1.0 


1.5 


1.85 


±0.13 


+0.15 
-0.11 


+0.08 
-0.07 




1.5 


2.0 


0.871 


±0.093 


+0.087 
-0.082 


+0.035 
-0.043 




2.0 


2.4 


0.74 


±0.11 


+0.07 
-0.15 


+0.04 
-0.03 


6 < < 9 


-1.5 


-1.0 


0.546 


±0.066 


+0.066 
-0.111 


+0.062 
-0.046 




-1.0 


-0.5 


1.374 


±0.081 


+0.166 
-0.043 


+0.097 
-0.065 




-0.5 


0.0 


2.067 


±0.097 


+0.200 
-0.163 


+0.098 
-0.091 




0.0 


0.5 


1.944 


±0.096 


+0.217 
-0.094 


+0.065 
-0.051 




0.5 


1.0 


1.429 


±0.093 


+0.184 
-0.059 


+0.048 
-0.040 




1.0 


1.5 


1.30 


±0.10 


+0.12 
-0.12 


+0.05 
-0.04 




1.5 


2.0 


0.542 


±0.075 


+0.071 
-0.080 


+0.009 
-0.027 




2.0 


2.4 


0.463 


±0.089 


+0.051 
-0.131 


+0.023 
-0.011 


4 ct > 9 


-1.50 


-1.00 












-1.0 


-0.5 


0.377 


±0.051 


+0.043 
-0.054 


+0.040 
-0.039 




-0.5 


0.0 


0.727 


±0.070 


+0.076 
-0.049 


+0.067 
-0.054 




0.0 


0.5 


0.950 


±0.074 


+0.130 
-0.045 


+0.072 
-0.070 




0.5 


1.0 


0.872 


±0.080 


+0.132 
-0.066 


+0.065 
-0.050 




1.0 


1.5 


0.559 


±0.079 


+0.087 
-0.071 


+0.032 
-0.028 




1.5 


2.0 


0.326 


±0.056 


+0.064 
-0.043 


+0.027 
-0.016 




2.0 


2.4 


0.276 


±0.063 


+0.049 
-0.066 


+0.013 
-0.019 



Table 2: The cross section da/drf Qt for events containing at least one D* meson 
in different regions of E 3 ^ 1 . The statistical (stat.), systematic (syst.) and energy 
scale (E-scale) uncertainties are shown separately. 



XXIII 



y t range 



da/dE^ ± stat. ± syst. ± El-scale (nb/GeV) 



_D*-tagged jet 



untagged jets 



6.00 9.00 

9.00 13.00 

13.00 18.00 

18.00 25.00 



0.889 
0.208 
0.0374 
0.0145 



±0.028 
±0.012 
±0.0059 
±0.0034 



+0.135 


+0.039 


-0.039 


-0.028 


+0.033 


+0.015 


-0.013 


-0.013 


+0.0064 


+0.0049 


-0.0054 


-0.0035 


+0.0046 


+0.0027 


-0.0043 


-0.0022 



0.704 
0.161 
0.0423 
0.0124 



±0.028 
±0.014 
±0.0076 
±0.0030 



+0.030 
-0.045 
+0.006 
-0.006 
+0.0122 
-0.0048 
+0.0040 
-0.0035 



+0.035 
-0.035 
+0.010 
-0.010 
+0.0049 
-0.0043 
+0.0007 
-0.0003 



Table 3: The cross section da/dE^ 1 for D*-tagged jets and untagged jets for 
—1.5 < rf ct < 2.4. The statistical (stat.), systematic (syst.) and energy scale 
(E-scale) uncertainties are shown separately. 



XXIV 



rf et range 



da/drf et ± stat. ± syst. ± E-scale (nb) 



D* matched jet 



other jets 



-1.50 

-1.00 

-0.50 

0.00 

0.50 

1.00 

1.50 

2.00 



-1.00 
-0.50 
0.00 
0.50 
1.00 
1.50 
2.00 
2.40 



0.412 
1.170 
1.831 
1.78 
1.18 
0.99 
0.19 



±0.055 
±0.075 
±0.097 
±0.11 
±0.11 
±0.13 
±0.17 



+0.060 
-0.054 
+0.236 
-0.064 
+0.324 
-0.123 
+0.34 
-0.12 

+0.30 
-0.08 
+0.24 
-0.11 
+0.22 
-0.17 



-0.043 
-0.034 
-0.089 
-0.048 
-0.084 
-0.067 
+0.05 
-0.04 

+0.03 
-0.03 
+0.02 
-0.01 
+0.00 
-0.02 



0.173 
0.560 
1.010 
1.113 
1.08 
1.014 
0.75 
0.70 



±0.043 
±0.075 
±0.094 
±0.095 
±0.10 
±0.098 
±0.10 
±0.13 



+0.058 
-0.097 
+0.039 
-0.029 
+0.081 
-0.224 
+0.099 
-0.056 
+0.10 
-0.07 
+0.085 
-0.129 
+0.10 
-0.11 
+0.08 
-0.20 



+0.022 
-0.015 
+0.035 
-0.036 
+0.053 
-0.063 
+0.044 
-0.044 

+0.05 
-0.04 
+0.054 
-0.046 
+0.01 
-0.04 
+0.04 
-0.02 



-1.50 
-1.00 

-0.50 

0.00 

0.50 

1.00 

1.50 

2.00 



-1.00 
-0.50 
0.00 
0.50 
1.00 
1.50 
2.00 
2.40 



0.386 
0.935 
1.367 
1.224 
0.771 
0.651 
0.113 



±0.052 
±0.060 
±0.073 
±0.074 
±0.069 
±0.086 
±0.100 



+0.051 
-0.051 
+0.157 
-0.037 
+0.186 
-0.060 
+0.182 
-0.075 
+0.162 
-0.050 
+0.132 
-0.076 
+0.154 
-0.098 



+0.040 
-0.032 
+0.071 
-0.038 
+0.063 
-0.050 
+0.037 
-0.024 
+0.019 
-0.017 
+0.015 
-0.010 
+0.001 
-0.013 



0.160 
0.438 
0.714 
0.725 
0.665 
0.623 
0.473 
0.466 



±0.040 
±0.058 
±0.067 
±0.062 
±0.063 
±0.061 
±0.066 
±0.089 



+0.055 
-0.089 
+0.039 
-0.023 
+0.055 
-0.148 
+0.065 
-0.036 
+0.055 
-0.039 
+0.049 
-0.080 
+0.064 
-0.070 
+0.049 
-0.131 



+0.021 
-0.014 
+0.027 
-0.028 
+0.037 
-0.045 
+0.029 
-0.028 
+0.033 
-0.025 
+0.033 
-0.028 
+0.008 
-0.022 
+0.023 
-0.011 



-1.50 
-1.00 

-0.50 

0.00 

0.50 

1.00 

1.50 

2.00 



-1.00 
-0.50 
0.00 
0.50 
1.00 
1.50 
2.00 
2.40 



0.239 ±0 
0.514 ±0 
0.620 
0.520 
0.231 
0.140 



±0 
±0 
±0 
±0 



036 
050 
055 
056 
057 
080 



+0.039 
-0.023 
+0.064 
-0.047 
+0.096 
-0.045 
+0.118 
-0.042 
+0.044 
-0.041 

+0.080 
-0.078 



+0.023 
-0.023 
+0.047 
-0.038 
+0.051 
-0.048 
+0.041 
-0.026 
+0.015 
-0.011 
+0.022 
-0.012 



0.148 
0.194 
0.328 
0.358 
0.335 
0.263 
0.275 



±0.043 
±0.055 
±0.051 
±0.063 
±0.055 
±0.048 
±0.062 



+0.021 
-0.091 
+0.058 
-0.015 
+0.053 
-0.018 
+0.038 
-0.052 

+0.067 
-0.039 
+0.048 
-0.033 
+0.049 
-0.065 



+0.021 
-0.019 
+0.018 
-0.015 
+0.020 
-0.021 
+0.023 
-0.026 
+0.016 
-0.018 
+0.021 
-0.013 
+0.013 
-0.019 



Table 4: The cross section da/drf Qt for D*-tagged jets and untagged jets for 
E 3 ^ 1 > 6 GeV. The statistical (stat.), systematic (syst.) and energy scale (E-scale) 
uncertainties are shown separately. 



XXV 



Pt(D*) range 


Ej^ range 


da/dE^ ± stat. ± syst. ± E-scale (nb/GeV) 


3.25 5.00 


6.0 9.0 
9.0 13.0 
13.0 18.0 
18.0 25.0 


o.907 ±0.032 I**" +g : g» 

0.137 ±0.015 $™ ±g : gJS 

0.0279 ±0.0085 j™« j™g 
0.0095 ±0.0047 ±glgggS ^oT 


5.00 8.00 


6.0 9.0 
9.0 13.0 
13.0 18.0 
18.0 25.0 


nein +nnifi +0-043 +0.004 
U.OIU ITU.U1U -0.039 -0.000 

1432 ±0 0091 +0 - 0127 +0 ' 0142 
u.i^oz zcu.uuyi _o.oo50 -0.0112 

0.0200 ±0.0040 ±8:88t? ±» 

0.0093 ±0.0034 Iffl Iffi* 


8.00 20.00 


6.0 9.0 
9.0 13.0 
13.0 18.0 
18.0 25.0 


n ns77 ±0 36 +0 0021 +-0 0018 

U.UOM JZU.UUOU -0.0071 0.0014 

o.o73i ±0.0050 ±8:881? tt 0009 

0.0338 ±0.0036 ^ ±8lS8S 
0.0092 ±0.0019 t° Zl + -° Zl 



Table 5: The cross section da/dEj, in bins of p®* The statistical (stat.), sys- 
tematic (syst.) and energy scale (E-scale) uncertainties are shown separately. 



XXVI 





ran ffp 


'/ 


ran f?p 


drr 1 dn }et 


± fit at 


± syst. 


-4- T^-sralp (nil] 


3.25 


5.00 


-1.5 


-1.0 


0.360 


±0.058 


+0.074 
-0.117 


+0.071 
-0.051 






-1.0 


-0.5 


0.903 


±0.075 


+0.282 
-0.105 


+0.113 
-0.083 






-0.5 


0.0 


1.465 


±0.096 


+0.162 
-0.184 


+0.140 
—0.122 






0.0 


0.5 


1.453 


±0.094 


+0.207 
-0.167 


+0.110 
-0.103 






0.5 


1.0 


1.063 


±0.097 


+0.244 
-0.046 


+0.094 
-0.074 






1.0 


1.5 


0.98 


±0.11 


+0.16 
-0.20 


+0.08 
-0.06 






1 .o 


9 n 




-+-n D7fi 


+0.078 
-0.048 


+0.021 
-0.025 






2.0 


2.4 


0.335 


±0.092 


+0.102 
-0.129 


+0.018 
-0.013 


5.00 


8.00 


-1.5 


-1.0 


0.211 


±0.033 


+0.081 
-0.051 


+0.014 
-0.007 






-1 


-0 5 


602 


±0 046 


+0.050 
-0.030 


+0.029 
-0.018 






-0.5 


0.0 


0.981 


±0.056 


+0.066 
-0.048 


+0.030 
-0.024 






0.0 


0.5 


1.015 


±0.057 


+0.105 
-0.074 


+0.026 
-0.022 






0.5 


1.0 


0.745 


±0.051 


+0.107 
-0.037 


+0.022 
-0.016 






1.0 


1.5 


0.534 


±0.053 


+0.052 
—0.050 


+0.018 
—0.014 






1.5 


2.0 


0.231 


±0.033 


+0.033 
-0.018 


+0.011 
-0.009 






2.0 


2.4 


0.240 


±0.045 


+0.035 
-0.049 


+0.011 
-0.009 


8.0 


20.0 


-1.50 


-1.00 














-1.0 


-0.5 


0.101 


±0.020 


+0.009 
-0.021 


+0.001 
-0.002 






-0.5 


0.0 


0.217 


±0.024 


+0.010 
-0.017 


+0.002 
-0.002 






0.0 


0.5 


0.274 


±0.027 


+0.019 
-0.011 


+0.002 
-0.001 






0.5 


1.0 


0.296 


±0.028 


+0.015 
-0.036 


+0.003 
-0.001 






1.0 


1.5 


0.178 


±0.024 


+0.008 
-0.018 


+0.003 
-0.003 






1.5 


2.0 


0.108 


±0.021 


+0.013 
-0.020 


+0.000 
-0.004 






2.0 


2.4 


0.092 


±0.020 


+0.014 
-0.019 


+0.004 
-0.005 



Table 6: The cross section da/drf ot in bins of plj! ■ The statistical (stat.), sys- 
tematic (syst.) and energy scale (E-scale) uncertainties are shown separately. 



XXVII 



x 7 (D*,jet) 


range 


da/dXj(D*, jet) 


± stat. 


± syst. 


± E-scale (nb) 


0.000 


0.250 


1.23 


±0.24 


+0.13 
-0.36 


+0.03 
-0.06 


0.250 


0.375 


2.15 


±0.28 


+0.37 
-0.30 


+0.09 
-0.11 


0.375 


0.500 


2.85 


±0.31 


+0.36 
-0.25 


+0.09 
-0.10 


0.500 


0.625 


3.82 


±0.30 


+0.46 
-0.36 


+0.04 
-0.06 


0.625 


0.750 


5.71 


±0.32 


+0.23 
-0.57 


+0.01 
-0.09 


0.750 


1.000 


2.82 


±0.15 


+0.25 
-0.25 


+0.37 
-0.30 



Table 7: The cross section da /cfe° bs (-D*, other jet). The statistical (stat.), sys- 
tematic (syst.) and energy scale (E-scale) uncertainties are shown separately. 



™obs 
x-y 


range 


da/dx° hs 


± stat. 


± syst. 


± E-scale (nb) 


0.000 


0.250 


0.19 


±0.15 


+0.34 
-0.07 


+0.02 
-0.01 


0.250 


0.375 


0.78 


±0.22 


+0.26 
-0.23 


+0.08 
-0.08 


0.375 


0.500 


0.99 


±0.19 


+0.24 
-0.23 


+0.08 
-0.09 


0.500 


0.625 


1.52 


±0.22 


+0.14 
-0.33 


+0.12 
-0.12 


0.625 


0.750 


2.03 


±0.24 


+0.21 
-0.23 


+0.18 
-0.16 


0.750 


1.000 


4.35 


±0.21 


+0.41 
-0.16 


+0.28 
-0.26 



Table 8: The dijet cross section da/dx^, for events containing at least one 
D* meson. The statistical (stat.), systematic (syst.) and energy scale (E-scale) 
uncertainties are shown separately. 



XXVIII 



A0J j 


range 


da/dAfti 


± stat. 


± syst. ± E-scale 


(nb/rad.) 


< x° bs < 1 


0.00 


1.57 


0.020 


±0.011 


+0.012 
-0.012 


+0.001 
-0.003 




1.57 


2.09 


0.149 


±0.035 


+0.051 
-0.046 


+0.025 
-0.014 




2.09 


2.36 


0.318 


±0.074 


+0.057 
-0.092 


+0.028 
-0.029 




2.36 


2.62 


0.67 


±0.10 


+0.17 
-0.10 


+0.07 
-0.05 




2.62 


2.88 


1.50 


±0.14 


+0.18 
-0.12 


+0.11 
-0.11 




2.88 


3.14 


4.06 


±0.22 


+0.39 
-0.32 


+0.24 
-0.23 




x° bs > 0.75 


0.00 


1.57 


0.0158 


±0.0068 


+0.0058 
-0.0115 


+0.0029 
-0.0017 




1.57 


2.09 


0.036 


±0.014 


+0.006 
-0.028 


+0.007 
-0.003 




2.09 


2.36 


0.126 


±0.038 


+0.032 
-0.027 


+0.011 
-0.012 




2.36 


2.62 


0.389 


±0.055 


+0.077 
-0.021 


+0.036 
-0.028 




2.62 


2.88 


0.829 


±0.091 


+0.105 
-0.066 


+0.055 
-0.055 




2.88 


3.14 


2.68 


±0.16 


+0.28 
-0.13 


+0.14 
-0.13 




if < 0.75 


0.00 


1.57 


0.0066 


±0.0093 


+0.0141 
-0.0074 


+0.0004 
-0.0010 




1.57 


2.09 


0.115 


±0.034 


+0.049 
-0.029 


+0.017 
-0.012 




2.09 


2.36 


0.190 


±0.062 


+0.042 
-0.082 


+0.018 
-0.018 




2.36 


2.62 


0.280 


±0.083 


+0.157 
-0.104 


+0.030 
-0.020 




2.62 


2.88 


0.67 


±0.11 


+0.10 
-0.10 


+0.06 
-0.06 




2.88 


3.14 


1.37 


±0.15 


+0.14 
-0.26 


+0.11 
-0.10 





Table 9: The dijet cross section dojd<$\ for events containing at least one D* 
meson for all x° bs , and for direct- enriched (i* s > 0.75,) and resolved- enriched 
^obs < 0.75,) samples. The statistical (stat.), systematic (syst.) and energy scale 
(E-scale) uncertainties are shown separately. 



XXIX 



Pj- range 


da/dpi^ ± stat. ± syst. ± E-scale (nb/GeV) 


< x° hs < 1 


3 
3 6 
6 9 
9 12 
12 18 


0.282 ±0.018 185S 
0.188 ±0.014 ±S;°» ^ 
0.0762 ±0.0095 ™& 18582 
0.0351 ±0.0079 1852? ±g:gggg 

o.oio9 ±0.0036 ig^ 18:888? 


if > 0.75 


3 
3 6 
6 9 
9 12 
12 18 


o.203 ±0.013 ±8:8fl ±8:8il 

0.1074 ±0.0092 t° S ±88811 

o.osoo ±0.0049 185888 ±8:88?i 
o.oi4i ±0.0045 ±8:881? 185858 

0.0066 ±0.0018 18582 ct8o°o°3 03 


x° hs < 0.75 


3 
3 6 
6 9 
9 12 
12 18 


o.o78 ±0.011 ±§;8?1 ±8:881 
o.o82 ±0.011 18:8?? 18:88? 
0.0479 ±0.0087 1™ ±8:881 
0.0207 ±0.0063 18:881 ±S 
o.oo46 ±0.0031 I858i ±8:8882 



Table 10: The dijet cross section da/dp^, for events containing at least one 
D* meson for all x° bs , and for direct- enriched fx° bs > 0.75,) and resolved- enriched 
(x° bs < 0.75,) samples. The statistical (stat.), systematic (syst.) and energy scale 
(E-scale) uncertainties are shown separately. 



XXX 



m JJ range 



da/dmP ± stat. ± syst. ± E-scale (nb/GeV) 



< x° bs < 1 



15 


20 


0.1299 


±0.0088 


+0.0095 
-0.0124 


+0.0080 
-0.0076 


20 


25 


0.0823 


±0.0076 


+0.0082 
-0.0076 


+0.0060 
-0.0059 


25 


30 


0.0484 


±0.0062 


+0.0106 
-0.0052 


+0.0049 
-0.0023 


30 


35 


0.0328 


±0.0056 


+0.0048 
-0.0080 


+0.0025 
-0.0043 


35 


40 


0.0223 


±0.0048 


+0.0018 
-0.0054 


+0.0033 
0.0009 


40 


50 


0.0081 


±0.0022 


+0.0017 
-0.0016 


+0.0006 
-0.0010 



x° bs > 0.75 



15 


20 


0.0788 


±0.0056 


+0.0064 
-0.0065 


+0.0033 
-0.0035 


20 


25 


0.0470 


±0.0049 


+0.0061 
-0.0029 


+0.0042 
-0.0028 


25 


30 


0.0319 


±0.0045 


+0.0049 
-0.0024 


+0.0026 
-0.0015 


30 


35 


0.0221 


±0.0045 


+0.0042 
-0.0045 


+0.0018 
-0.0031 


35 


40 


0.0136 


±0.0035 


+0.0019 
-0.0045 


+0.0022 
0.0003 


40 


50 


0.0040 


±0.0016 


+0.0021 
-0.0013 


+0.0003 
-0.0004 



x° bs < 0.75 



15 


20 


0.0512 


±0.0069 


+0.0055 
-0.0097 


+0.0051 
-0.0043 


20 


25 


0.0357 


±0.0060 


+0.0043 
-0.0083 


+0.0016 
-0.0034 


25 


30 


0.0166 


±0.0041 


+0.0059 
-0.0036 


+0.0023 
-0.0008 


30 


35 


0.0109 


±0.0033 


+0.0030 
-0.0043 


+0.0007 
-0.0013 


35 


40 


0.0086 


±0.0033 


+0.0022 
-0.0027 


+0.0011 
0.0006 


40 


50 


0.0038 


±0.0015 


+0.0009 
-0.0013 


+0.0004 
-0.0006 



Table 11: The dijet cross section da/dM^, for events containing at least one 
D* meson for all a;° bs , and for direct- enriched (x^ s > 0.75) and resolved- enriched 
^obs < 0.75,) samples. The statistical (stat.), systematic (syst.) and energy scale 
(E-scale) uncertainties are shown separately. 
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ZEUS 




0.14 0.15 0.16 0.17 

AM (GeV) 

Figure 1: The distribution of the mass difference, AM = M(Ktttt s ) — M{Kir), 
for D* candidates with a single jet. The D*^ candidates (dots) are shown compared 
to the wrong charge combinations (histogram). The dashed vertical lines show the 
signal region, 0.143 < AM < O.lJf.8 GeV. The number of D* mesons is determined 
by subtracting the wrong charge background as described in Section IJ.il The fit is 
for illustrative purposes only. 
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Figure 2: The distribution of AR between D* mesons and each jet in the event. 
The data (dots) are compared to Herwig (solid line) and Pythia (dashed line) 
MC predictions. The MC predictions are area normalised to the data using the 
normalisation factors shown in the figure. The dotted vertical line indicates the 
AR = 0.6 cut which separates D* -tagged jets from untagged jets. The AR(D*,jet) 
distribution for D*-tagged jets for the MC hadron-level predictions is shown for 
Herwig (thin solid line) and Pythia (thin dashed line). The hadron-level predic- 
tions are area normalised to the data in the region AR(D*,jet) < 0.6 
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Figure 3: Cross-section da /dE^ for jets in events (dots) containing at least one 
D* meson for different regions in rj^ ct . The comparison is made to massive QCD 
predictions with (solid line) and without (dashed line) hadronisation corrections 
applied. The theoretical uncertainties (hatched band) come from the change in 
scales simultaneously with the change in charm mass. The beauty component is 
also shown (lower histogram) . 
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Figure 4: Cross-section da/dE^ 1 for D* -tagged jets and untagged jets (dots). 
The comparison is made to massive QCD predictions with (solid line) and without 
(dotted line) hadronisation corrections applied. The beauty component is also shown 
(lower histogram) . For the untagged jet distribution, the massless QCD predictions 
are also shown with (solid line) and without (dashed line) hadronisation corrections 
applied. The theoretical uncertainties (hatched bands) come, in the case of the 
massive calculations, from changing renormalisation and factorisation scales as 
well as the charm mass simultaneously. In the case of the massless calculations, 
they come from changing the scales only. The prediction with no component from 
b-quark fragmentation to a D* is also shown (dotted line). 
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Figure 5: Cross-section da/drf et for jets in events (dots) containing at least one 
D* meson for different regions in . The comparison is made to massive QCD 
predictions with (solid line) and without (dotted line) hadronisation corrections 
applied. The theoretical uncertainties (hatched band) come from the change in 
scales simultaneously with the change in charm mass. The beauty component is 
also shown (lower histogram) . 
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Figure 6: Cross-section da/drf ct for D* -tagged jets and untagged jets (dots). 
The comparison is made to massive QCD predictions with (solid line) and with- 
out (dashed line) hadronisation corrections applied. The beauty component is also 
shown (lower histogram) . For the untagged jet distribution, the massless QCD pre- 
dictions are also shown with (solid line) and without (dashed line) hadronisation 
corrections applied. The theoretical uncertainties (hatched bands) come, in the case 
of the massive calculations, from changing renormalisation and factorisation scales 
as well as the charm mass simultaneously. In the case of the massless calculations, 
they come from changing the scales only. The prediction with no component from 
b-quark fragmentation to a D* is also shown (dotted line). 
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Figure 7: Cross-sections da / dE^ and da / drf et in bins of p®* ■ The data (solid 
dots) are compared to the massive QCD predictions with (solid line) and with- 
out (dotted line) hadronisation corrections applied. The theoretical uncertainties 
(hatched band) come from the change in scales simultaneously with the change in 
charm mass. The beauty component is also shown (lower histogram) . 
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theoretical uncertainties (hatched bands ) come, in the case of the massive calcula- 
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Figure 9: Dijet cross-sections (a) da/dx° bs , (b) da/dA^, (c) da/d(p^) 2 and (d) 
da/dM^ for the data (solid dots) compared to massive QCD predictions with (solid 
line) and without (dotted line) hadronisation corrections applied. The theoretical 
uncertainties (hatched band) come from the change in scales simultaneously with 
the change in charm mass. The beauty component is also shown (lower histogram) . 



XL 



> o.i 
o 

>< 

+ 

T 0.05 

* 

Q 
+ 

T 



73 





0.1 



ZEUS 

o.i 



x° bs >0.75 



a 



ZEUS 98-00 

Jet energy scale uncertainty 



20 



30 



40 



50 



x° bs <0.75 



0.05 



NLO QCD (massive) 

NLO QCD (massive) ® had. 

Beauty 



0.05 







-i 1 1 1 r- 



-i 1 1 1 1 r- 



; x° bs >0.75 



C 



HERWIGx 2.5 
PYTHIAx 1.5 



U 



_i I i i i_ 



j i I i i i_ 



20 30 40 50 

M jj (GeV) 




0.05 



40 50 

M jj (GeV) 



obs 



> 



Figure 10: Cross-section da/dM^ separated into (a,c) direct enriched (x\ 
0.75,) and (b,d) resolved enriched fx° bs < 0.75,). The data (solid dots) are compared 
(a,b) to the massive QCD prediction with (solid line) and without (dotted line) 
hadronisation corrections applied. The theoretical uncertainties (hatched band) 
come from the change in scales simultaneously with the change in charm mass. 
The beauty component is also shown (lower histogram). The data are also com- 
pared (c,d) with Herwig (solid line) and Pythia (dashed line) MC predictions 
multiplied by the indicated factors. 
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Figure 11: Cross-section da/dAcf) 33 separated into (a,c) direct enriched (x° > 
0.75 and (b,d) resolved enriched fx° bs < 0.75). The data (solid dots) are compared 
(a,b) to the massive QCD prediction with (solid line) and without (dotted line) 
hadronisation corrections applied. The theoretical uncertainties (hatched band) 
come from the change in scales simultaneously with the change in charm mass. 
The beauty component is also shown (lower histogram). The data are also com- 
pared (c,d) with Herwig (solid line) and Pythia (dashed line) MC predictions 
multiplied by the indicated factors. 
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Figure 12: Cross-section da/d(p^) 2 separated into (a,c) direct enriched (x° bs > 
0.75) and (b,d) resolved enriched fx° bs < 0.75). The data (solid dots) are compared 
(a,b) to the massive QCD prediction with (solid line) and without (dotted line) 
hadronisation corrections applied. The theoretical uncertainties (hatched band) 
come from the change in scales simultaneously with the change in charm mass. 
The beauty component is also shown (lower histogram). The data are also com- 
pared (c,d) with Herwig (solid line) and Pythia (dashed line) MC predictions 
multiplied by the indicated factors. 
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